function run_ball(animalID)
[xls_num,xls_txt]=xlsread('\\argon.fmi.ch\keller.g\ExpLog\Stacks.xlsx');
col = find(strcmp(xls_txt(1,:), 'AnimalID') == 1);
animal_inds=find(strcmp(xls_txt(2:end,col),animalID));
% % % stck_col=find(strcmp(xls_txt(1,:), 'StackID') == 1);
exp_col=find(strcmp(xls_txt(1,:), 'ExpID') == 1);
% % % StckID=xls_num(animal_inds,stck_col);
ExpID=unique(xls_num(animal_inds,exp_col));

cutExp=input(['Cut Exp (1 or 0)?']);
if cutExp
    cutExp_ind=input(['Which ones [enter as vector]: ']);
    for tnd=1:length(cutExp_ind)
        ExpID(ExpID==cutExp_ind(tnd))=[];
    end
end

auxdata={};
maxX=[];
sess_mat=[];
for fnd=1:length(ExpID)
    curr_exp=ExpID(fnd);
    curr_stck=xls_num(xls_num(:,2)==curr_exp,3);
    curr_data_all=[];
    for dnd=1:length(curr_stck)
        curr_data=load_lvd(['\\KELLER-RIG1-ANA\RawData\zmarpawe' '\' animalID '\' 'S1-T' num2str(curr_stck(dnd)) '.lvd']);
        display(['Loading ' num2str(curr_stck(dnd))]);
        curr_data=curr_data(5,:);
        if dnd==1
            strt=0;
        else
            strt=length(curr_data_all);
        end
        curr_data_all(1,strt+1:(strt+length(curr_data)))=curr_data(1,:);
        sess_mat(fnd,dnd)=strt+length(curr_data);
    end
    auxdata{fnd}=curr_data_all;
    maxX(fnd)=length(auxdata{fnd});
end
max_x=max(maxX);
% % % figure;
% % % for znd=1:length(ExpID)
% % %     hold on
% % %     subplot(length(ExpID),1,znd)
% % %     plot(auxdata{znd})
% % %     xlim([0 max_x]);
% % %     title(num2str(ExpID(znd)))
% % % end
vel_mat={};
perc=[];
for gnd=1:length(ExpID)
    vel=diff(auxdata{gnd});
    vel(vel>5)=vel(vel>5)-10;
    vel(vel<-5)=vel(vel<-5)+10;
    velM_raw=vel;
    vel=ftfil(vel,1000,0,10);
    vel_ind=((smooth2(abs(vel),200))>0.01);
    vel_mat{gnd}=vel_ind;
    perc(gnd)=(sum(vel_ind)/length(vel_ind))*100;
end

figure;
for znd=1:length(ExpID)
    hold on
    subplot(length(ExpID),1,znd)
    plot(vel_mat{znd})
    set(gca,'xlim',([0 max_x]),'XTick',[0:1000000:max_x],'ylim',([0 2]),'YTick',([0:1:2]))
    title([num2str(ExpID(znd)) ' ' num2str(round(perc(znd))) '% time spent running'])
    for hnd=1:length(sess_mat(1,:))
        if sess_mat(znd,hnd)~=0
            line([sess_mat(znd,hnd) sess_mat(znd,hnd)],[0 2],'LineStyle',':','color','r')
        end
    end
end
